


// IV Regressions by time periods. 1-5 6-10-11-16


////////////////////////////////////////////////
//   time=[1-5]                               //
////////////////////////////////////////////////



use "G:\Data\Workdata\707677\Common\01Data\02STATA\panel_iv", clear 

cd  G:\Data\Workdata\707677\Common\03Output\02Regresssions





keep if time==1

tab obs_lawyer

egen lawyer_id_unique =group(lawyerid court) // Unique lawyer id within courts. some lawyers work in multiple courts
bys lawyer_id_unique: egen n_lawyer=count(lawyer_id_unique)
egen mean_lawyer_leniency= mean(granted), by(lawyer_id_unique)

gen iv_judge2_all = (mean_lawyer_leniency*n_lawyer - granted)/(n_lawyer-1)

gen iv_lawyer2=iv_judge2_all-(iv_court-granted)/(obs_court-1)

* Data restrictions more than 20 obs pr lawyer and multiple lawyers in a given court*year



keep if n_lawyer>=20

egen n_judge=nvals(lawyer_id_unique), by(court aar)


tab n_judge

keep if n_judge>1

keep pnr iv_lawyer2 lawyer_id_unique


merge 1:m pnr using "G:\Data\Workdata\707677\Common\01Data\02STATA\panel_iv"


keep if _merge==3
drop _merge


keep if time>=1 & time<=5

tab court, gen(courtdum)

tab time, gen(timedum)


egen court_pyear =group(court pyear) 
tab court_pyear

tab aar, gen(aardum)
drop aardum1



local Xvar_short "male age_0_40 age_41_50 age_51_60 age_61_70 erhvervsindk_mean employed_pstill_mean unemployed_pstill_mean legmarried_m1 racedum1 housedum1 passiv_m1 aktiv_m1 higher_secondary university social_mean edu_missing housing_missing timelon_dum1 timelon_dum2 timelon_dum3 timelon_dum4 house_m1"


local depvar "employed_pstill unemployed_pstill out_of_lf_pstill timelon_real aktiv_real passiv_real formue_real koejd_real house tilbtot_real disability kontant_real social netovskud_real dispon_ny_real legmarried attrition mortgage_real bankgaeld_real"



tab time

/* IV Regress */


reg   erhvervsindk_real granted `Xvar_short' i.court_pyear if time==1

local individuals=e(N)	
 

sum erhvervsindk_real if granted==0

local dismissed_mean=r(mean) 

ivreghdfe erhvervsindk_real (granted=iv_lawyer2) `Xvar_short' aardum*, absorb(court_pyear) tol(1e-6) first cluster(lawyer_id_unique) 


outreg2 using iv_results_1_5.xls, replace ctitle(IV erhvervsindk_real) keep(granted) ///
 addstat(Non-granted mean, `dismissed_mean',individuals, `individuals')



foreach var of varlist `depvar'{
	
reg   `var' granted `Xvar_short' i.court_pyear if time==1
local individuals=e(N)	

sum `var' if granted==0
local dismissed_mean=r(mean)

ivreghdfe `var' (granted=iv_lawyer2) `Xvar_short' , absorb(court_pyear) tol(1e-6) first cluster(lawyer_id_unique) 


outreg2 using iv_results_1_5.xls, append ctitle(IV `var')  


}

////////////////////////////////////////////////
//   time=[6-10]                               //
////////////////////////////////////////////////



use "G:\Data\Workdata\707677\Common\01Data\02STATA\panel_iv", clear 

cd  G:\Data\Workdata\707677\Common\03Output\02Regresssions



keep if time==1

tab obs_lawyer

egen lawyer_id_unique =group(lawyerid court) // Unique lawyer id within courts. some lawyers work in multiple courts
bys lawyer_id_unique: egen n_lawyer=count(lawyer_id_unique)
egen mean_lawyer_leniency= mean(granted), by(lawyer_id_unique)

gen iv_judge2_all = (mean_lawyer_leniency*n_lawyer - granted)/(n_lawyer-1)

gen iv_lawyer2=iv_judge2_all-(iv_court-granted)/(obs_court-1)


* Data restrictions more than 20 obs pr lawyer and multiple lawyers in a given court*year

keep if n_lawyer>=20

egen n_judge=nvals(lawyer_id_unique), by(court aar)


tab n_judge

keep if n_judge>1

keep pnr iv_lawyer2 lawyer_id_unique
/

merge 1:m pnr using "G:\Data\Workdata\707677\Common\01Data\02STATA\panel_iv"


keep if _merge==3
drop _merge


keep if time>=6 & time<=10


tab court, gen(courtdum)

tab time, gen(timedum)


egen court_pyear =group(court pyear) 
tab court_pyear

tab aar, gen(aardum)
drop aardum1 



local Xvar_short "male age_0_40 age_41_50 age_51_60 age_61_70 erhvervsindk_mean employed_pstill_mean unemployed_pstill_mean legmarried_m1 racedum1 housedum1 passiv_m1 aktiv_m1 higher_secondary university social_mean edu_missing housing_missing timelon_dum1 timelon_dum2 timelon_dum3 timelon_dum4 house_m1"


local depvar "employed_pstill unemployed_pstill out_of_lf_pstill timelon_real aktiv_real passiv_real formue_real koejd_real house tilbtot_real disability kontant_real social netovskud_real dispon_ny_real legmarried attrition mortgage_real bankgaeld_real"



tab time

/* IV Regress */


reg   erhvervsindk_real granted `Xvar_short' i.court_pyear if time==6


local individuals=e(N)	
 
 
sum erhvervsindk_real if granted==0

local dismissed_mean=r(mean) 


ivreghdfe erhvervsindk_real (granted=iv_lawyer2) `Xvar_short' aardum*, absorb(court_pyear) tol(1e-6) first cluster(lawyer_id_unique) 


outreg2 using iv_results_6_10.xls, replace ctitle(IV erhvervsindk_real) keep(granted) ///
addstat(Non-granted mean, `dismissed_mean',individuals, `individuals')



foreach var of varlist `depvar'{
	
reg   `var' granted `Xvar_short' i.court_pyear if time==6
local individuals=e(N)	

sum `var' if granted==0
local dismissed_mean=r(mean)

ivreghdfe `var' (granted=iv_lawyer2) `Xvar_short' , absorb(court_pyear) tol(1e-6) first cluster(lawyer_id_unique) 


outreg2 using iv_results_6_10.xls, append ctitle(IV `var')  




}


////////////////////////////////////////////////
//   time=[11-16]                             //
////////////////////////////////////////////////



use "G:\Data\Workdata\707677\Common\01Data\02STATA\panel_iv", clear 

cd  G:\Data\Workdata\707677\Common\03Output\02Regresssions



keep if time==1

tab obs_lawyer

egen lawyer_id_unique =group(lawyerid court) // Unique lawyer id within courts. some lawyers work in multiple courts
bys lawyer_id_unique: egen n_lawyer=count(lawyer_id_unique)
egen mean_lawyer_leniency= mean(granted), by(lawyer_id_unique)

gen iv_judge2_all = (mean_lawyer_leniency*n_lawyer - granted)/(n_lawyer-1)

gen iv_lawyer2=iv_judge2_all-(iv_court-granted)/(obs_court-1)

* Data restrictions more than 20 obs pr lawyer and multiple lawyers in a given court*year


keep if n_lawyer>=20

egen n_judge=nvals(lawyer_id_unique), by(court aar)


tab n_judge

keep if n_judge>1

keep pnr iv_lawyer2 lawyer_id_unique



merge 1:m pnr using "G:\Data\Workdata\707677\Common\01Data\02STATA\panel_iv"


keep if _merge==3
drop _merge

keep if time>=11 & time<=16

tab court, gen(courtdum)

tab time, gen(timedum)


egen court_pyear =group(court pyear) 
tab court_pyear

tab aar, gen(aardum)
drop aardum1




local Xvar_short "male age_0_40 age_41_50 age_51_60 age_61_70 erhvervsindk_mean employed_pstill_mean unemployed_pstill_mean legmarried_m1 racedum1 housedum1 passiv_m1 aktiv_m1 higher_secondary university social_mean edu_missing housing_missing timelon_dum1 timelon_dum2 timelon_dum3 timelon_dum4 house_m1"



local depvar "employed_pstill unemployed_pstill out_of_lf_pstill timelon_real aktiv_real passiv_real formue_real koejd_real house tilbtot_real disability kontant_real social netovskud_real dispon_ny_real legmarried attrition mortgage_real bankgaeld_real"




tab time

/* IV Regress */


reg   erhvervsindk_real granted `Xvar_short' i.court_pyear if time==11
local individuals=e(N)	
 

sum erhvervsindk_real if granted==0

local dismissed_mean=r(mean) 


ivreghdfe erhvervsindk_real (granted=iv_lawyer2) `Xvar_short' aardum*, absorb(court_pyear) tol(1e-6) first cluster(lawyer_id_unique) 


outreg2 using iv_results_11_16.xls, replace ctitle(IV erhvervsindk_real) keep(granted) ///
addstat(Non-granted mean, `dismissed_mean',individuals, `individuals')



foreach var of varlist `depvar'{
	
reg   `var' granted `Xvar_short' i.court_pyear if time==11
local individuals=e(N)	

sum `var' if granted==0
local dismissed_mean=r(mean)

ivreghdfe `var' (granted=iv_lawyer2) `Xvar_short', absorb(court_pyear) tol(1e-6) first cluster(lawyer_id_unique) 


outreg2 using iv_results_11_16.xls, append ctitle(IV `var')  




}


